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Abstract 

Within the framework of quantum electrodynamics (QED), vacuum is a nonhnear medium which 
can be hnearized for a rapidly time-varying electromagnetic field with a small amplitude subjected to 
a magnetostatic field. The linearized QED vacuum is a uniaxial dielectric-magnetic medium for which 
the degree of anisotropy is exceedingly small. By implementing an affine transformation of the spatial 
coordinates, the degree of anisotropy may become sufficiently large as to be readily perceivable. The 
inverse Bruggeman formalism can be implemented to specify a homogenized composite material (IfCM) 
which is electromagnetically equivalent to the affinely transformed QED vacuum. This HCM can arise 
from remarkably simple component materials; for example, two isotropic dielectric materials and two 
isotropic magnetic materials, randomly distributed as oriented spheroidal particles. 
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1 Introduction 

Classical vacuum is a linear medium for which the principle of superposition holds. Consequently, light 
propagation in classical vacuum is unaffected by the presence of a magnetostatic field. However, within the 
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framework of quantum electrodynamics (QED), vacuum is a nonlinear medium [1]. The QED vacuum can 
be linearized for a rapidly time- varying electromagnetic field with a small amplitude subjected to a slowly 
varying (or static) magnetic field [2] . A consequence of linearization is that the QED vacuum appears as a 
uniaxial dielectric-magnetic medium for optical fields [3]. The constitutive parameters which characterize 
this uniaxial medium depend on the magnitude and direction of the magnetostatic field. 

The degree of anisotropy associated with the QED vacuum is exceedingly small. Consequently, a direct 
measurement of this attribute poses enormous challenges to experimentalists [4] , and experimental verifica- 
tion of the anisotropy of the QED vacuum is eagerly awaited [5]. In view of this difficulty, we propose an 
experimental simulation of the QED vacuum which would enable the anisotropy to be explored for prac- 
ticable magnetostatic fields. The simulation is based on a homogenized composite material (HCM), which 
arises from the homogenization of remarkably simple component materials. For example, the component 
materials could be isotropic dielectric and magnetic materials, randomly distributed as oriented spheroidal 
particles. Similar HCM-based simulations have recently been described for the Schwarzschild-(anti-)de Sitter 
spacetime [6] and cosmic strings [7]. 

In the following sections, 3-vectors are underlined with the addition of the'symbol denoting a unit vector. 
Double underlining indicates a 3x3 dyadic with the dyadic transpose being labelled with the additional 
symbol. The 3x3 identity dyadic represented as /. The permittivity and permeability of classical vacuum 
are written as eo = 8.854 x 10~^^ F m~^ and /xo = 47r x 10""^ H m~^, respectively. 

2 Electrodynamics of QED vacuum 

We consider vacuum under the influence of a magnetostatic field S^. = In classical vacuum, the 

passage of light is unaffected by B^, as reported by an inertial observer. However, this is not the case 
for the QED vacuum. The QED vacuum is a nonlinear medium which can be linearized for rapidly time- 
varying plane waves. Thereby, for propagation of light, QED vacuum is represented by the anisotropic 
dielectric-magnetic constitutive relations [3] 



The relative permittivity and relative permeability dyadics of the QED vacuum have the uniaxial forms [8] 



where the constant a = 6.623 x 10^^^ H^^ kg^^ m^ s^. The constitutive dyadics (2) were derived by Adler 
[3] from the Heisenbcrg-Eulcr effective Lagrangrian of the electromagnetic field [9, 10]. 

Since a is exceedingly small, the value of a\B_^\'^ is also exceedingly small in comparison with unity 
for typical values of |^^|. For example, in the PVLAS experiment \B_J = 5 T typically [5], which yields 
= 1.656 X 10~^^. Accordingly, the degree of anisotropy represented by the constitutive dyadics (2) 
is also exceedingly small. In order to achieve degrees of anisotropy that could be realistically attained in a 
controlled manner for a practical simulation of QED vacuum, we implement the affine transformation 




(1) 




X x' — J • X 



(3) 



of the spatial coordinates. The transformation dyadic 




(4) 
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employs 



9 = 



/(l-8a|Bj2)(i + 20a|Sj2) 
Y {1-2S\BJ){1 + 55\BJ) 

l-8a|Sj2 



(5) 



and the scalar parameter 6 > 0. For definiteness, we fix i5 = 0.02. Thus, the affine-transformed relative 
permittivity and permeability dyadics are given as [11] 



1 



J '6' .r 



det e 

(l-25|BJ)(l-44 



and 



1 



det /i — = — 

1-2S\BJ 



2 V= ^s±=Ls 



(l-8a|BJ2) 



1 - 4a|Sj2 (1 + 120a|Bj2) ---^ 



(6) 

(7) 
(8) 

(9) 

(10) 
(11) 



Notice that for the range \B_J G (0, 10) T, the denominators of both terms on the right side of Eq. (10) are 
both approximately equal to unity, and therefore e' « /i'. The components e'^ and are linearly dependent 

upon I-BJ, as illustrated in Fig. 1. The plots of /i'^ ^ versus jS^j are practically identical to those of e'^ 



3 Simulation as a homogenized composite material 

Let us now turn to the question: How can one specify an HCM which is a uniaxial dielectric-magnetic 
material with relative permittivity dyadic ^^jqj^j = d and relative permeability dyadic /i = /^'? In 
order to answer this question, we make use of the well-established Bruggeman homogenization formalism 
[12, 13]. 

Suppose we consider the homogenization of four component materials, labelled a, 6, c and d. Two of 
the components (a and b, say) are isotropic dielectric materials while the other two (c and d) arc isotropic 
magnetic materials. Thus, the component materials are specified by 

(i) the relative permittivities ea, Ec and e^, with Ec = = 1; and 

(ii) the relative permeabilities /io, A*f), /^c and /id, with jia = fit = ^■ 

The four component materials are randomly distributed with volume fractions fa, ft, fa fd G (0, 1), with 
fd = ^ — fa — fb ^ fc- AH four component materials consist of identically oriented spheroidal particles. The 
symmetry axis for all these spheroidal particles lies parallel to B^^ . Accordingly, the surface of each spheroid 
relative to its centre is prescribed by the vector 

rs=Pig^-L, (12) 

wherein the shape dyadic 

= (l - 4s,s) + UeB^B^, {£ = a, b, c, d) , (13) 
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is real symmetric [14] and positive definite, the radial unit vector is f , and the linear measure pi is required 
to be small compared to the electromagnetic wavelengths under consideration. The shape parameter J/^ > 1 
for prolate spheroids and Q <Ui < \ for oblate ones. 

Under the Bruggeman homogenization formalism, the corresponding HCM is a uniaxial dielectric-magnetic 
material, specified by relative permittivity and permeability dyadics of the form 



= HCM ~ ^* 



(l - lA) + r^^'^ksB,, (r = (14) 



For the particular case of the uniaxial dielectric-magnetic HCM involved here, full details of the numerical 
process of computing the dyadics e and ^ , from a knowledge of ea.b,c,d, fJ-a,b.c,d, Ua.b,c,d and fa,b.c.d, 
are provided elsewhere [6] . 

Conventionally, homogenization formalisms are used to estimate the constitutive parameters of HCMs, 
based on a knowledge of the constitutive and morphological parameters of their component materials and 
their volume fractions. In contrast, here our goal is to estimate the constitutive and morphological parameters 
as well as the volume fractions of the component materials which give rise to a HCM such that e ^^j^^ coincides 
with e' and n coincides with u'. We do so via an inverse implementation of the Bruggeman formalism. 

= =HCM = 

Formal expressions of the inverse Bruggeman formalism are available [15], but in some instances these can be 
ill-defined [16]. In practice, the inverse formalism may be more effectively implemented by direct numerical 
methods [17]. Note that certain constitutive parameter regimes have been identified as problematic for the 
inverse Bruggeman formalism [18], but these regimes are not the same as those considered here. 

We consider the following three different implementations of the inverse Bruggeman formalism. In each 
implementation, four scalar parameters are to be determined. 

I. The relative permittivities ea.t and the relative permeabilities fic,d sltc assumed to be known, and all 
spheroidal particles have the same shape, i.e., Ua ^ Ub = Uc ~ Ud ^ U. Wc then determine the 
common shape parameter U and the volume fractions /„, ft and fe- 
ll. The relative permittivities ea,& and the relative permeabilities fic,d sue assumed to be known, and the 
volume fractions fa.b,c arc fixed. We then determine the shape parameters Ua, Ub, Uc and Ud- 

111. The shape parameters Ua,b,c,d and the volume fractions fa,b.c are fixed. We then determine the relative 
permittivities Ca.b and relative permeabilities /ic,rf. 

To describe the inversion of the Bruggeman formalism, let us focus on implementation I as a repre- 
sentative example, the inversion processes for implementations II and III being analogous. Suppose that 
{e^'-^^"'^, ef^'^*'', /if'^'-^*^, /i^*-^^'^} are forward Bruggeman estimates of the HCM's relative permittivity and 
relative permeability parameters which are computed for physically reasonable ranges of the parameters U 

and /a,&,c; i.e., U G (C/",C/+) and fa.b,c e [fa,b.cJtb.^- 

(1) Let fa = if a + /+) /2, fb = {fb + fb) /2, and = (/- + /+) /2. For aU U e (C/^, C/+), determine 
the value which yields the minimum value of the scalar quantity 



s \ I / till t s Ms \ , I t^t n 



1/2 

■ (15) 



(2) Let U = U\ fb = {fb + fb) /2, and = (/- + /+) /2. For all e (/-,/+), determine the value /] 
which yields the minimum value of A. 

(3) Let U = C/t, /, = fl and /, = (/- +/+)/2. For ah fb G {f^J^)^ determine the value fl which 
yields the minimum value of A. 

(4) Let U = W , fa = fl, and fb = fl. For all fc £ (fcifc)' determine the value /] which yields the 
minimum value of A. 
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The steps (1)~(4) are then repeated, with fl, fl, and fl being the fixed values of /a,6,c in step (1), fl and 
fl being the fixed values of f^^c in step (2), and fl being the fixed value of /c in step (3), until the value of 
A becomes acceptably small. 

4 Numerical illustrations 

Numerical illustrations of the implementations I-III are provided in Figs. 2-4. For all results presented, the 
degree of convergence of the numerical schemes which provide the inverse Bruggeman estimates was < 1%, 
and in most instances this value was < 0.1%. 

I. For Fig. 2, the constitutive parameters of the component materials were taken to be = 4, ef, = 0.3, 
He = 3.4 and = 0.4. The computed common shape parameter U and volume fractions fa,b,c are 
plotted versus \B_^\. While the volume fractions vary little as \B_J\ is increased from 1 to 2.5 T, the 
common shape parameter increases exponentially. 

II. The constitutive parameters of the component materials were again taken to be = 4, Cb = 0.3, 
He = 3.4 and = 0.4 for Fig. 3. In addition, the volume fractions were fixed at fa = 0.15, fb — 0.25 
and fc — 0.21. The computed four shape parameters Ua.b,c,d are plotted versus \B_J. All four shape 
parameters increase uniformly as \B_J is increased from 1 to 2 T. This refiects the fact that the degree 
of anisotropy of the HCM is required to increase as \B_J increases. 

III. Lastly, in Fig. 4 the common shape parameter is fixed at U = 5 while the volume fractions arc fixed at 
fa.b,c = 2.5. The computed constitutive parameters ea,& and fj,c,d are plotted versus l-B,]. In this case, 
Ca turns out to be approximately the same as ^c- And similarly turns out to be approximately the 
same as /i^- While and /ic decrease uniformly as \B_g \ is increased from 1 to 3 T, the opposite is true 
of Cb and Hd ■ 

5 Closing remarks 

By means of the inverse Bruggeman formalism, an HCM may be specified which is electromagnetically 
equivalent to the QED vacuum subject to a spatial affinc transformation. The affinely transformed QED 
vacuum retains the same uniaxial dielectric-magnetic form as the un-transformcd QED vacuum, but the dcgce 
of anisotropy is greatly exaggerated by means of the affinc transformation. By reversing the transformation 
represented by eq. (3), the properties of QED vacuum may be inferred from those of the HCM. 

For illustration, the inverse homogenization formulation presented here was based on four isotropic com- 
ponent materials. However, the desired HCM could also be realized by alternative inverse homogenization 
formulations. For example, the HCM could arise from only two component materials. These two compo- 
nents materials could be either both isotropic dielectric-magnetic materials distributed as oriented spheroidal 
particles or both uniaxial dielectric-magnetic materials (with parallel symmetry axes) distributed as spher- 
ical particles [13, 19]. However, the four-component formulation presented here involves the simplest of 
component materials and allows a large degree of freedom in choosing their constitutive parameters. 

Finally, let us note that the relative permittivities and relative permeabilities of the component materials 
needed for the HCM, as presented in Figs. 2-4, are not at all infeasible. Indeed, present-day technology allows 
for the possibility of materials with a considerably wider range of constitutive parameters to be engineered 
[20, 21, 22]. 

Acknowledgment: AL thanks the Charles Godfrey Binder Endowment at Pcnn State for partial financial 
support of his research activities. 
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I I (T) 

Figure 1: The relative permittivity parameters (blue, dashed) and (red. solid) plotted versus |_B J (T). 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

I ^ l(T) 

Figure 2: Implementation I. The common shape parameter U (thick solid, red) and volume fractions fa 
(dashed, green), fb (broken dashed, blue), and /c (thin solid, blue) plotted versus \B_g\ (T). The relative 
permittivities = 4, = 0.3, /ic = 3.4 and = 0.4. 
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1.0 1.2 1.4 1.6 1.8 2.0 

I ^ l(T) 



Figure 3: Implementation II. The shape parameters Ua (thick solid, red), Ut (dashed, green), Uc (broken 
dashed, blue), and Ud (thin solid, blue) plotted versus \B_^\ (T). The relative permittivities ~ 4, e;, = 0.3, 
fic = 3.4 and fid — 0.4; and the volume fractions fa ~ 0.15, fb = 0.25 and fc = 0.21. 




I ^ l(T) 

Figure 4: Implementation III. The relative permittivities (thick solid, red) and ef, (dashed, green) and the 
relative permeabilities /.ic (broken dashed, blue) and fid (thin solid, blue) plotted versus \B_g\ (T). The shape 
parameter U = Ua,b,c,d = 5 and volume fractions fa,b,c = 0.25. 
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